(N 



^D 



Oh 
CD 



o 

(N 



c^ 



IISER/HEP/03/12 

Many faces of low mass neutralino dark matter in the 
unconstrained MSSM, LHC data and new signals 

Arghya Choudhury 1 and Amitava Datta 2 

Indian Institute of Science Education and Research - Kolkata, 

Mohanpur Campus, PO: BCKV Campus Main Office, 



pg \ Nadia, West Bengal - 741252, India. 



Abstract 

If all strongly interacting sparticles (the squarks and the gluinos) in an uncon- 



t^i \ strained minimal super symmetric standard model (MSSM) are heavier than the corre- 

p+ 

sponding mass lower limits in the minimal supergravity (mSUGRA) model, obtained 



by the current LHC experiments, then the existing data allow a variety of electroweak 
(EW) sectors with light sparticles yielding dark matter (DM) relic density allowed 

(n: 

£> . by the WMAP data. Some of the sparticles may lie just above the existing lower 

CO ' bounds from LEP and lead to many novel DM producing mechanisms not common 



in mSUGRA. This is illustrated by revisiting the above squark-gluino mass limits ob- 
tained by the ATLAS Collaboration, with an unconstrained EW sector with masses 
not correlated with the strong sector. Using their selection criteria and the correspond- 
ing cross section limits, we find at the generator level using Pythia, that the changes 
rN ■ in the mass limits, if any, are by at most 10-12 % in most scenarios. In some cases, 



however, the relaxation of the gluino mass limits are larger (~ 20%). If a subset of the 
strongly interacting sparticles in an unconstrained MSSM are within the reach of the 
LHC, then signals sensitive to the EW sector may be obtained. This is illustrated by 
simulating the blj fop, I = e and u, and brjpT signals in i) the light stop scenario and 
ii) the light stop-gluino scenario with various light EW sectors allowed by the WMAP 
data. Some of the more general models may be realized with non-universal scalar and 
gaugino masses. 
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1 Introduction 

The ATLAS and the CMS collaborations have been searching for supersymmetry (SUSY) 
[1] at the ongoing LHC 7 TeV experiments [2, 3, 4, 5]. No signal has been seen so far. A 
number of phenomenological analyses of the prospect of susy search at 7 Tev have also been 
published [6]. 

The recently announced results for 1 fb~ x data [2, 3, 4, 5] have been presented as con- 
straints in the popular minimal supergravity (mSUGRA) model [7] 3 . As expected the jets 
+ fix channel, which arises from squark-gluino pair production in all combinations with 
large cross sections and is not suppressed due to small branching ratios, yields the strongest 
constraints on squark-gluino masses. 

However, due to the special correlations among the superparticle (sparticle) masses in 
mSUGRA, the above mass bounds on strongly interacting sparticles impose stringent indirect 
lower bounds on the sparticle masses in the electroweak (EW) sector consisting of the sleptons 
and the electroweak gauginos. In many cases these model dependent bounds are significantly 
stronger than the corresponding direct bounds from LEP [8]. On the other hand, as we 
shall elaborate below by revising the existing limits, for any choice of squark-gluino masses 
compatible with the jets + fix data there may exist a variety of EW sectors with much 
lighter sparticles compared to that in mSUGRA. In fact the jets + fix data is only mildly 
sensitive to the electroweak sector in most cases (some exceptions will be listed below). Thus 
it is desirable to think of signals at the LHC which are directly sensitive to the electroweak 
sector. 

One of the main attractive features of R-parity conserving SUSY is that the lightest 
supersymmetric particle (LSP) is stable. In many models the weakly interacting lightest 
neutralino (xi) * s assumed to be the LSP and it turns out to be a very popular candidate 
for the observed dark matter (DM) in the universe [9, 10]. The observed dark matter (DM) 
relic density (Qh 2 ) in the universe has been precisely measured by the Wilkinson Microwave 
Anisotropy Probe (WMAP) collaboration [11] and if 10% theoretical uncertainty is added 
[12] then DM relic density is bounded by 0.09 < Vth 2 < 0.13 at 2a level. EW sectors with 
relatively low mass sparticles can provide many DM producing mechanisms yielding relic 
densities consistent with the data. As mentioned in the last paragraph such possibilities are 



3 The model we consider is also referred to as the constrained minimal supersymmetric standard model 
(CMSSM) in the literature. 



now mostly excluded in the mSUGRA framework by the LHC experiments. The negative 
impact of early LHC data on low mass neutralino DM and on the prospect of direct DM 
search experiments were noted in [13]. However, this exclusion based on strongly model 
dependent assumptions is certainly not the final verdict on neutralino DM since the SUSY 
breaking mechanism is essentially unknown. Hence it is worthwhile to revisit the viable DM 
producing mechanisms within frameworks more general than mSUGRA invoking as little 
model dependent assumptions as is practicable. 

It has recently been emphasized [14] that there are many relic density producing mech- 
anisms involving light electroweak sectors, which are practically independent of the strong 
sector leaving aside the possibility of the coannihilation of the LSP with the lighter top 
squark (t\) [15]. Thus any model with all strongly interacting sparticles beyond the reach of 
the LHC - 7 TeV experiments and a relatively light electroweak sector consistent with the 
observed relic density is allowed by the LHC data. In [14] no model dependent correlation 
among the sparticle masses in the strong and EW sectors was imposed. However, for the 
sake of simplicity it was assumed that the masses in the EW sector are correlated as in 
mSUGRA. Several observable signals involving parameter spaces consistent with both LHC 
and WMAP data were proposed. In this paper we give up the last assumption and study 
the impact of light EW sectors in an unconstrained minimal supersymmetric extension of 
the standard model (MSSM) on both LHC signals and the relic density data. 

Obviously the best way to test such models would be signals at a new e + — e~ collider 
directly sensitive to only EW sparticles with masses ~ the corresponding lower limits from 
LEP. In the absence of such an accelerator one can look,e.g., for the clean 31 (1 = e or /i) 
signal from Chargino -neutralino pair production [16] and di-lepton + $t signal from slepton 
pair production [17]. Apriority both the signals are viable even if the strongly interacting 
sparticles are heavy. 

The experience from the simulations of the LHC-14 TeV experiments [18], however, 
does not encourage optimistic expectation at the on going experiments. In the former case 
chargino-neutralino masses modestly above the LEP limits were found to be observable. One 
can improve the chargino-neutralino mass reach to some extent by considering 21 + 1 r and 
11 + 2 r events along with the 31 events [19]. It is also estimated that sleptons with masses 
S 400 GeV can be probed at the LHC-14 TeV with C £ 30 fb^ 1 (see [18]). 

Models in which only a sub-set of the strongly interacting sparticles along with light EW 



sparticles are within the reach of the current LHC experiment are also compatible with the 
data. The decay of the strongly interacting sparticles into final states involving EW sparticles 
may provide new signals of reasonable size. Some examples illustrating this [3, 14, 20] have 
already been discussed. However, the analyses in [14, 20] were done when either the 1 
fb _1 data were not available or available in the unpublished form without many details. In 
this paper we shall improve these analyses using the published data and more general EW 
sectors consistent with the WMAP data as elaborated above. Our main goals are to check 
the impact of the more general models on i) the squark-gluino mass limits already obtained 
within the framework of mSUGRA and ii) the viability of the signals proposed in [14, 20] 
for more general scenarios. 

In the simplest scenario considered in this paper following [14, 20], the lighter stop 
squark t\ is assumed to be the only strongly interacting sparticle accessible to the current 
LHC experiments. Other groups have also investigated the light t\ scenarios [22, 23]. If 
this squark is the next lightest supersymmetric particle (NLSP), then it will decay via the 
loop induced final state consisting of a charm quark and the LSP. As noted earlier the DM 
relic density may be produced in such a scenario via ti-LSP coannihilation. This channel 
has been investigated recently but does not appear to be very promising [24, 25]. Recently 
novel signals of t± NLSP have been proposed [26, 27]. However, the competition between 
the above mode and the four body decay of the t± NLSP [28] may further complicate the 
issue. The consequences of this competition in the context of the Tevatron experiments have 
already been discussed [29]. 

It was recently pointed out in ref [14, 20] that if this squark decays dominantly into 
the lighter chargino (x^) an d a b-quark, then viable signals sensitive to the EW sector may 
appear in the blj fix and brj fix channels. However, for ti within the reach of the ongoing 
LHC experiments the above signals cannot survive the strong cuts on $p and M e ff usually 
employed for general squark-gluino searches in jets + $t by the LHC collaborations. Softer 
dedicated cuts need to be employed for this signal as shown in [14, 20]. In this paper we 
shall revisit the signals for more general EW sectors as elaborated above. 

Another class of models interesting in the context of the LHC experiments are the ones 
with both the t\ and the gluinos are within the reach of the current experiments: the 
light stop-gluino (LSG) model [14, 23]. This scenario can be realized in models with non- 
universal scalar and gaugino masses at the GUT scale [14]. The same signals as discussed 



in the last paragraph may arise in this case also via the production of gluino pairs which 
then decay dominantly into tti pairs. It may be recalled that the ATLAS collaboration has 
already investigated the blj fix signal [3]. However, the hard cuts proposed for separating 
the background eliminates the events from the direct t\ pair production. As a result the 
gluino mass limit (500-520 GeV) obtained by them is by and large independent of the vti-r . 
On the other hand the alternative selection criteria with softer cuts as proposed in [14, 20] 
can separate the blj fix events from t\ pair production from the ones coming from gluino 
pair production. In this paper we shall follow this approach and shall revisit the signal for 
the more general EW sectors discussed above. Moreover the complementary channel brj ft? 
which is more important when the electroweak gauginos decay dominantly in r-rich final 
states was not considered in [3]. This signal was studied for the first time in [14]. 

The plan of the paper is as follows. Section 2 contains all the results obtained in this 
paper and the related discussions. In Section 2.1 we list the input parameters for the uncon- 
strained EW sector. The main departures from the mSUGRA spectrum are duly emphasized. 
The possibility of accommodating these departures with non-universal scalar and gaugino 
masses are discussed as and when appropriate. In sub-section 2.2 the parameter space in 
the unconstrained EW sector consistent with WMAP data is delineated. DM relic density 
producing mechanisms which are not allowed in mSUGRA but occur frequently in more 
general scenario are pointed out. Revision of the squark-gluino mass bounds in mSUGRA 
for more general EW sectors consistent with WMAP data is presented in the sub-section 
2.3. Novel signal in the blj Pt and brj fix channels in the light stop and LSG scenarios are 
also included in this section. Our conclusions will be summarized in Section 3. 



2 Results and Discussions : 

2.1 Models and Parameters : 

The popular mSUGRA model [7] has only five free parameters including soft SUSY breaking 
terms. These are mo (the common scalar mass), Tn\i% (the common gaugino mass), A$ (the 
common trilinear coupling), all given at the gauge coupling unification scale {Mq ~ 2 x 10 16 
GeV); the ratio of the Higgs vacuum expectation values at the electroweak scale namely tan/3 
and the sign of //. The magnitude of \x is determined by the radiative electroweak symmetry 



breaking condition. 

In contrast in the phenomenological MSSM Mi,M2, M3 are the three gaugino mass 
parameters at the weak scale and no special relation among them is assumed. It is known 
for a long time that even in mGUGRA type models where masses at the weak scale are 
determined by some boundary conditions at a high scale, a variety of gaugino mass relations 
at the weak scale may emerge due to non-universal boundary conditions for gaugino masses 
at M G [30]. 

We also consider the following weak scale parameters. Mq. is the mass parameter for 
the ith (i = 1, 2 and 3 is the generation index) generation squarks of the L-type belonging 
to a doublet of SU(2)l- Ms. (Mj.) is the mass parameter for the R-type singlet up (down) 
squarks. Similarly M~ v (M? t ) is the left (right) type slepton mass parameter for the ith 
generation. A t , A^, A T are the 3rd generation trilinear soft couplings. The other input 
parameters are the mass of pseudoscalar Higgs boson rriA, the higgsino mixing parameter \x 
and tan/3. 

In our numerical analysis we have assumed that all squarks and gluinos are beyond the 
reach of the LHC 7 TeV run 4 . For the sake of simplicity we have taken M^ = Mq 2 = Mq 3 
= M^ = Ma 2 = M^ 3 = Mj = Mj = Mj = 1.5 TeV. For the relic density computation 
in this subsection we also take M3 (which defines the gluino mass) = 1.5 TeV at the weak 
scale. 

We have further assumed that Mj = Mf = M fl = M f2 = Mj ; Mj = M? 3 ; but these 
common masses at the weak scale are treated as free parameters. The equality of the slepton 
mass parameters at the weak scale is not realized in mSUGRA, where L-type sleptons are 
typically heavier than the R-type. Since the sneutrino mass is correlated with that of L-type 
sleptons, it also exceed R- type slepton masses over most of the parameter space. On the 
other hand models with non-universal scalar masses at the high scale can yield considerably 
different mass hierarchies among the two types of sleptons at the weak scale. The running 
of the common scalar masses at the SUSY breaking scale (say, the Planck scale) to the 
GUT scale (Mq) may create the above nonuniversality when L and R type sleptons belong 
to different representations of the GUT group [31]. When a GUT group breaks down to a 
group of lower rank, certain U(l) symmetry breaking D-terms can also lead to such non- 
universality [32]. 



'All masses and parameters having the dimension of mass arc in GeV unless stated otherwise. 



Relatively light sleptons and, consequently, sneutrinos can affect the DM relic density 
production through LSP annihilation via slepton exchange and / or via LSP-sneutrino co- 
annihilation in a large region of the parameter space where the sneutrino happens to be 
the next lightest supersymmetric particle (NLSP). In the context of collider phenomenology, 
the relative probabilities of m-lepton + n-jet + fix signatures, for different m and n, may 
turn out to be quite distinct in comparison to mSUGRA [33] in such models. This may 
distinguish between different models. 

In addition we assume A t = A^ = A T = -600; /x = 362.0 ; m^ = 1 TeV and tan/3 = 
10. Moreover since our focus is on a light EW sector the mass parameters in this sector are 
restricted - somewhat arbitrarily - to be less than 200. As a result we have either a bino 
dominated LSP or an admixture of bino and wino with negligible Higgsino component. 

In this parameter space the computed mass of the lighter Higgs scalar (m^), after taking 
into account a theoretical uncertainty of 3 [34], satisfies the LEP bound m^ > 114.4 [35]. 
Further discussion of rrih in the light of more recent experiments is given below. 

We have used micrOMEGAs (v.2.4.1) [36] for computing the DM relic density. The susy 
particles spectra and the decay branching ratios (BRs) have been computed by SUSPECT 
[37] and SDECAY [38]. 

In Fig. 1 we plot Qh 2 vs Mi for different choices of the other parameters in the EW 
sector. We have also checked that varying the mj or mg does not affect the DM relic density, 
unless the difference between mj and m~o is so small that the possibility of stop-neutralino 
coannihilation [15] opens up. 

To begin with we have searched for the parameter space allowed by the WMAP data. 
Results are given in Fig. 1. Here we have plotted Qh 2 vs Mi for different M2 and common 
slepton masses. We have chosen the following five representative sets of M 2 and the common 
slepton mass of the third generation (Mj- = M? 3 ): model-1 (110,115), model-2 (110,150) , 
model-3 (150,115), model-4 (150,200) and model-5 (200,200). The choices M 2 = 110, 150 
and 200 yield m~± = 106 ,145 and 191 respectively while Mj 3 = M H = 115, 150 and 200 
result in m~ = 89, 131, 186 respectively. For each combination of M 2 and M ?3 , we have 
varied the common slepton mass (Mj) for the first two generations in steps of 10 in the range 
100 to 200. Thus we allow the lighter f mass eigenstate (t\) to be lighter or heavier than 
the sleptons belonging to the first two generations. This is also a significant departure from 
the mSUGRA spectrum. It should be noted that we have several choices where the masses 



of the sparticles in the electroweak sector are just above the corresponding LEP limits. 
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Figure 1: Mi vs flh 2 for different choices of M 2 , M H (= M~ h ) and M r (= M~ h = M u = M fl = M? 2 ). Blue 
horizontal lines represent the WMAP allowed DM relic density band (0.09 to 0.13). 

In Fig 1, each color or point shape corresponds to a particular choice of M2 and M? 3 as 
indicated in the upper right corner. The ten different lines of the same color or type in Fig. 
1 correspond to different choice of Mj. 

Throughout this work the pole mass of the top quark (running bottom quark mass 



evaluated in the MS scheme) will be taken as m t (ra&) = 173.2 (4.25). We have also 
included the current constraints on squarks, sleptons, gluino, chargino, neutralinos and the 
lighter scalar Higgs boson masses obtained from collider data. For example in the CMSSM 
models with tan/3 = 10, A = and \x > 0, ATLAS has excluded squarks and gluinos having 



equal mass below 950 GeV [2]. We shall assume that m~± > 103.5, all the sleptons (except 
the lighter stau mass eigenstate) are heavier than 100. This is basically a simplified form of 
the LEP limits. The lighter f mass eigenstate is assumed to be heavier than 86 as required 
by the LEP data. We emphasize that there is no model independent bound on m~o from 

X\ 

collider data. From the chargino mass bound from LEP there is a bound m~o > 50.8 for 

Xi 

tan/3 = 10 and // > 0, which is valid in mSUGRA. From LHC, WMAP and XENON100 [39] 
there is a bound m~o > 160 in mSUGRA [40]. We have only considered parameter spaces 
resulting in 114.4 [35] < m^ < 127 [41, 42]. We have checked that recent hint of a peak at 
nih ~ 125 reported by the LHC collaborations [41], which is not yet statistically significant, 
can be accommodated if we make Aq large (say 1.5 TeV). However, our parameter space 
allowed by the relic density data is not drastically affected by the choice of Aq. The impact 
of the rrih bound on different SUSY models and various observables has been considered by 
several groups [43]. 

2.2 Summary of Fig.l : 

There are five sets of data : 

• Model- 1 : The red points (marked by +) represent the computed Qh 2 in this param- 
eter space. Here both Xi and T\ lie just above the LEP limit. The latter is the NLSP 
for the entire parameter space scanned and its mass sets the limit for the LSP mass. 

• Model-2 : The squares represent the computed relic density. Here the mass of the Xi 
is just above the LEP limit and it is the NLSP for the whole parameter space studied. 
The mass of f\ is somewhat away from the LEP limit. 

• Model-3 : The green cross marks represent the parameter space with M 2 = 150, M? 3 
= 115. Here the mass of f\ is just above the LEP limit and it is the NLSP in our study. 
However, the mass of Xi is away from the LEP limit. 

• Model-4 : Relic densities are represented by the violet pentagons. In this case both 
Xi and T\ have masses far away from the the LEP limit. For relatively low common 
slepton masses the sneutrino is the NLSP. Otherwise the chargino is the NLSP. 



• Model-5 : In this case the orange triangles represent the relic densities. Over most of 
the parameter space the sneutrino is the NLSP. For the highest common slepton mass 
(200), the lighter stau, the sneutrino and the chargino have closely spaced masses. 

We have also demarcated the WMAP allowed DM relic density band (0.09 to 0.13) by 
the blue horizontal lines. In the next section we will take different benchmark points from 
this region for studying the impact of this more general EW sector on LHC physics. 

The importance of LSP annihilation via relatively low mass R-type slepton exchange 
(Bulk annihilation) in producing the observed relic density is well known in mSUGRA [9, 10]. 
The corresponding region of the parameter space (the Bulk region) is, however, strongly dis- 
favoured by the current LHC constraints as they do not permit the sleptons to be sufficiently 
light. In the unconstrained MSSM this mechanism retains it importance for certain ranges 
of the LSP mass where both L and R type slepton exchange may be important. However, 
even for smaller and larger LSP masses there are many other potentially important relic 
density producing mechanisms. 

In order to streamline the discussions, we divide the whole region allowed by the WMAP 
data into three parts according to the value of M\ as follows : 

I) 40 ~ Mi ~ 60 : In this region different choices of M 2 consistent with LEP data give 
almost the same DM relic density, which, however, depends strongly on the slepton masses. 
The neutralinos studied in this paper are either bino dominated or mixture of bino and wino. 
For such neutralinos, m~o significantly smaller than 40 are not allowed due to the LEP lower 
bounds on chargino and slepton masses and the WMAP constraints. In fact we have checked 
that m~a ~ 35 in the parameter space we have scanned. 

Three main mechanisms or combinations operate in this case: i) bulk annihilation, ii) 
annihilation mediated by a Z boson of low virtuality or iii) mediated by the lighter Higgs 
scalar (h) of low virtuality. It should be noted that even if we ignore the indirect bound on 
m~o from the LHC bounds in mSUGRA, a large part of this region is strongly disfavoured 
in mSUGRA due to the LSP mass bound from LEP inferred from the chargino mass limit. 
We also stress that no co-annihilation process can play any role in this parameter space due 
to the LEP lower bounds on EW sparticle masses. 

For stau mass a little above the LEP bound (Model 1 or 3), M- { not very far from the 
corresponding bound and the LSP mass below the Z resonance, this happens through pure 
bulk annihilation for a range of slepton masses(e.g., with M x = 40, Mj <130). For higher 
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Mi on both sides of the Z resonance, both bulk annihilation and annihilation via a virtual 
Z contributes adequately in all models. However, the slepton masses for which this happens 
is model dependent. For Mi on the Z resonance or in its immediate vicinity the annihilation 
cross section is too large and Vth 2 is too small for all LEP allowed slepton and stau masses. 

Well above the Z resonance, the Z contribution begins to reduce. Now process hi) takes 
over. Combination of i) and iii) is now the potential mechanism. This happens in all models 
1-5 but the range of Mj is again model dependent. There is even a small parameter space 
where i), ii) and iii) all contributes significantly. 

On the h resonance or in its immediate neighbourhood Qh 2 is again too small. At Mi 
= 56, Qh 2 falls abruptly and nearly at Mi = 59 we see a minima. This is due to the higgs 
resonance (here computed rrih = 113 ± 3 (theoretical uncertainty). Also here the relative 
contribution to l/(Qh 2 ) is largest from the process x?Xi -^ bb (nearly 80 %). 

II) 60 ~ Mi ~ 80 : A large number of points in model 1 and 3, with t\ NLSP, for different 
common slepton masses lie in the region allowed by WMAP data. The value of Vth 2 remains 
almost constant for LSP masses in this range for fixed values of other parameters. Pure bulk 
annihilation is the source of DM relic density. 

In models 2, 4 and 5 with heavier f\, pure bulk annihilation is not viable. However, if the 
sneutrinos belonging to the first two generations are the NLSP, then a small but significant 
contribution from LSP- sneutrino coannihilation along with bulk annihilation may produce 
the required relic density. For each Mi there is a narrow range of sneutrino masses where 
this happens. This mechanism is indicated by a sharp fall of the relic density over a narrow 
range of Mi. 

III) 80 ~ Mi ~ 200 : In models 1 and 3, the LSP-ri coannihilation becomes effective 
for a small range of Mi just above 80. 

In models 2, 4 and 5 with heavier fi, bulk annihilation alone can not produce the 
observed relic density. If the lighter chargino is the NLSP, several choices of the common 
slepton mass yield approximately the same relic density (see, e.g., the squares (model 2) 
in the neighbourhood of Mi ~ 95). In fact Qh 2 in this case primarily depends on the 
LSP and chargino properties alone. LSP annihilation into W pairs and LSP-Gaugino co- 
annihilation are the significant mechanisms for relic density production. If, on the other 
hand, the sneutrino is the NLSP, then sneutrino co-annihilation and bulk annihilation may 
serve the purpose. This happens for a small range of the common slepton mass for each Mi. 
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Several examples are shown in Fig. 1. 

2.3 Collider Signals 

In this section we discuss the collider signatures corresponding to some benchmark points 
allowed by both WMAP (see Section 2.2) and LHC data. 

In this paper all leading order (LO) signal cross-sections have been computed by CalcHEP 
[44] unless otherwise stated. For any two body final state (except for QCD processes) with 
identical particles or sparticles both the renormalization and the factorization scales are 
taken as, [ar = /af = M, where M is the mass of the particle or sparticle concerned 5 . 
For two unequal masses in the final state the scales are taken to be the average of the two. 
For QCD events the scales have been chosen to be equal to vl which is the energy in the 
parton CM frame, and the cross- section is computed by Pythia [45]. All LO cross- sections 
are computed using CTEQ5L parton density functions (PDFs) [46]. 

The next to leading order (NLO) cross-sections for the signal processes have been com- 
puted by PROSPINO [48] using the CTEQ5M PDFs. The K-factors are computed by 
comparing with the LO cross-section. The LO cross- sect ions from PROSPINO agree well 
with CalcHEP for the same choice of the scales. The NLO cross sections are used to compute 
revised mass limits and the signal sizes presented in this section. 

The NLO background cross-sections are not known for some backgrounds - in particu- 
lar for the QCD processes. For computing the significance of the signal we conservatively 
multiply the total LO background by an overall factor of two. 

We have considered the backgrounds from ti, QCD events , W + n-jets events and Z + 
n-jets, where W and Z decays into all channels. After the final cuts ti, W + lj and W + 2j 
are the main backgrounds, ti events are generated using Pythia and the LO cross-section 
has been taken from CalcHEP which is 85.5 pb. QCD processes are generated by Pythia in 
different p T bins : 25 < p T < 400, 400 < p T < 1000 and 1000 < p T < 2000 , where p T is 
defined in the rest frame of the parton collision. The main contribution comes from the low 
Pt bin, which has a cross-section of ~ 7.7E + 07 pb. However, for the other two bins the 
background events are negligible. 

For W + n-jets events we have generated events with n — 0,1 and 2 at the parton level 



5 See footnote 4. 
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using ALPGEN (v 2.13) [47]. We have generated these events subjected to the condition 
that Pj, > 20, AR(j,j) > 0.3 and \rj\ < 4.5. These partonic events have been fed to Pythia 
for parton showering, hadronization, fragmentation and decays etc. 

We have used the toy calorimeter simulation (PYCELL) provided in Pythia with the set- 
tings described in [14, 20]. Lepton (7 = e, fi) selection, b- jet tagging and r- jet identification 
are also implemented following these references. 

The strongest limits on squark, gluino masses come from the jets + fix channel. The 
ATLAS group has introduced five sets of selection criteria (SC) (see Table 1 of [2]) for new 
physics search in this channel. The observed number of events in this channel and the SM 
background estimated from the data for each of the above SC lead to a model independent 
upper bound on a new , the effective cross-sections for any new physics scenario. Five such 
bounds, thus obtained, are 22 fb, 25 fb, 429 fb, 27 fb and 17 fb respectively at 95% confidence 
level. These bounds can be used to obtain exclusion contours in any specific SUSY breaking 
model. 

Some examples of the limits (L1-L4) obtained by the ATLAS collaboration in the 
mSUGRA/CMSSM model (See Fig. 2 (right) of [2]) for tan/3 = 10, A = and // > 
are given below. 

Ll)Squarks and gluinos of equal mass are excluded for masses below 950 GeV. Here the 
average squark mass is considered. 

From the exclusion plot of ATLAS it also follows that m = 1.0 (1.5) TeV and vti\ii ~ 
300 (210) lies on the exclusion contour. It then follows that 

L2) For average squarks mass of 1.1 TeV gluino masses below 780 GeV are excluded. 

L3) For average squarks mass of 1.5 TeV gluino masses below 600 GeV are excluded. 

For the validation our simulation, we focus on the above exclusion contour. By definition 
&new = crAe, where a is the raw production cross-section in any new physics scenario, A is 
the acceptance and e is the efficiency [2]. We have chosen several points on this exclusion 
contour. For each point we compute A and e corresponding to all the five SC of ATLAS 
from our simulation. For each mSUGRA point a in the NLO is computed by PROSPINO. 
We find that each computed o new is reasonably close to at least one of the upper limits (see 
above) obtained by ATLAS. This reflects that our simulation in mSUGRA would lead to an 
exclusion contour pretty close to the one obtained by ATLAS. We follow this procedure to 
obtain new limits in other models. 
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We now turn our attention to more general models where the masses of the sparticles 
in the strong and EW sectors are assumed to be uncorrelated. We consider various EW 
sectors allowed by the WMAP data as illustrated in Fig. 1 of this paper and compute the 
squark-gluino mass limits from ATLAS bounds on a new . 



Points 


M x 


M 2 


Mi 


Mr 




m xf 


Tu- 


mr 


PI 


45 


150 


160 


200 


43 


145 


rn 


166 


P2 


128 


150 


160 


200 


121 


145 


186 


166 


P3 


43 


150 


200 


115 


40 


145 


89 


205 


P4 


81 


150 


200 


115 


77 


145 


89 


205 


P5 


45 


110 


150 


150 


42 


106 


131 


157 


P6 


99 


150 


120 


200 


93 


145 


186 


128 


P7 


94 


110 


200 


150 


88 


106 


131 


205 



Table 1: Parameters and mass spectra corresponding to different benchmark points taken from 
different regions of Fig. 1. 

For a better understanding of the revised limits, we also select some benchmark EW 
sectors from Fig.l (see Table 1). The corresponding BRs relevant for this discussion are 
presented in Table 2. Some important characteristics of the selected points are given below. 

• The point PI corresponds to a low m~o not allowed by the LEP data in mSUGRA. 

• For the point P2, the BR of the mode xt — > Xi^u l — e an d /J>, is large (78 %) 
compared to what is obtained over most of the mSUGRA parameter space. This is 
due to the fact that the decay Xi ~ *" XiW is not allowed and the common slepton mass 
is relatively low. 

• For P3 and P4, the combined BR of the modes Xi ~~ ** ^\ v t an d v t t is 100 %. 

• In examples P5 and P7, xt is the NLSP. But the mass difference (Am) between the 
Xi and Xi is large in P5 (65)and much smaller in P7 (19). The BR of the mode 
xt -> X°M is 48 (23) % for the point P5 (P7). 
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In P6, vt is the NLSP, Xi is lighter than T\ but heavier than II and vg. Here BR of 
Xi ~~ >" vd is 87% and xf ~^ Il v i is 13 %■ 



Decay Modes 


PI 


P2 


P3 


P4 


P5 


P6 


P7 


Xi -> XiW 


- 


17 


- 


- 


30 


- 


51 


-> X?^/ 


- 


78 


- 


- 


48 


- 


23 


-»■ x?rz/ r 


- 


5 


- 


- 


22 


- 


26 


-^TiM,. 


- 


- 


36 


36 


- 


- 


- 


^Xi^ 


100 


- 


- 


- 


- 


- 


- 


—7- z/ r r 


- 


- 


64 


64 


- 


- 


- 


-> z^Z 


- 


- 


- 


- 


- 


87 


- 


-^IlVi 


- 


- 


- 


- 


- 


13 


- 


n -> X?r 


75 


49 


100 


100 


78 


52 


53 


-►xSr 


9 


21 


- 


- 


8 


19 


19 


~^xt V T 


16 


30 


- 


- 


14 


29 


28 



Table 2: The BRs (%) of the dominant decay modes of \\ an d t\ for the benchmark points. 

In order to compare the new limits with the ones in mSUGRA, we fix the average squarks 
mass as in the examples LI, L2, L3 given above and vary m-x. We then compute the jets 
+ fix signature corresponding to EW sectors randomly chosen from the points in Fig. 1 
allowed by WMAP data. If a gluino mass yields cross sections smaller than the ATLAS 
upper bounds on a new for all selection criteria, we reduce it keeping average squark mass 
fixed. In this way we obtain a relaxed limit on m-x for a fixed average squark mass when the 
cross-section exceeds the upper bound on a new [2] for at least one selection criteria. 

For most of the points considered by us, the m-x limit changes at most by 10 to 12 % 
compared to the mSUGRA limit. For e.g., for the benchmark point PI the lower bounds on 
m-x in LI, L2, L3 are reduced by 50 to 70. As noted earlier this is due to relatively low m~o. 

The leptonic BR of the Xi I s large for P6. The lepton veto in the ATLAS selection 
criteria relaxes the limits on m-x in all cases by 70-80. This is reminiscent of the relaxation of 
the m~- m-x mass limits obtained from Tevatron data (see the 1st paper of [33]) for similar 
reasons. 
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On the other hand if BR(x x — > V T r, T\U T ) = 100 % then the limit on m-x is strengthened 
by 5-10 % due to the large BR of hadronic r decays. The relative changes in the limits noted 
by us are typically of the order of the uncertainties in the production cross sections due to 
the choice of QCD scale and/or the parton density function. This simple exercise illustrates 
that ATLAS limits on squarks-gluino masses are approximately valid for many different 
electroweak sectors quite different from that in mSUGRA. However, we also arrive at the 
important conclusion that slepton and / or chargino masses just above the corresponding 
LEP limits are very much allowed by the current LHC data. 



ATLAS CUTS [2] 
'High Mass' 


mSUGRA point with 

m = 2500, rai/2 = 185 

A = and tan(3=10 


Benchmark point PI 

with m-x = 550 


P T of leading jet > 130 


92% 


76% 


fi T > 130 


54% 


33% 


P T of 2nd jet > 80 


75% 


68% 


P T of 3rd jet > 80 


79% 


78% 


P T of 4th jet > 80 


57% 


61% 


Pr/M eff > 0.2 


47% 


32% 


M eff > 1100 


29% 


18% 



Table 3: Efficiency table for different cuts employed by ATLAS (see Table 1 of [2]) for deriving in 
mSUGRA the limit L4 (see text) and the corresponding efficiencies for the benchmark point PI. 

L4) For heavy squarks (say 2.5 or 3.0 TeV) ATLAS data has excluded gluinos with 
masses in the range 550 to 570 . 

Now keeping the average squarks mass as in L4 but choosing the EW sectors from Figure 
1, we have checked the change in the limit on gluino mass. Some interesting examples are 
given below. Models with EW sector as in PI significantly relax the limit on the gluino mass 
to 450. This is the largest relaxation of m-x limit we have noted in this paper. In mSUGRA 
the limit L4 is obtained by the last cut set ('High Mass') employed by the ATLAS group [2]. 
The efficiencies for these cuts are presented in Table 3. For the same average squark mass 
and m-x, the corresponding efficiencies for the point PI are displayed in Table 3 which are 
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found to be smaller in each case. Models with EW sectors as in P2 (P3) relaxes (strengthens) 
the gluino mass limit to 500 (620). 

We next consider the light stop scenario (i.e., t\ is the only strongly interacting sparticle 
within the reach of 7 TeV experiments) and examine the potential of the blj fix and brj ft? 
signals stemming from t\ pair production for various EW sectors chosen from Fig 1. For 
this analysis we have taken the third generation squark mass parameters (Mq 3 ,Mn 3 and Mj ) 
as free parameters. For t\ masses within the reach of the LHC experiments at 7 TeV, this 
scenario is not constrained at all by the jets + fix data obtained with hard cuts on kinematic 
variables like fix and m e ff. 



EW sector 
same as 


S (blj fi T ) 
(Cut Set 1 of [20]) 


S (brjP T ) 
(Cut Set 1 of [U]) 


PI 


426(4.8*) 


82(2.0**) 


P2 


449(5.0) 


8(0.2) 


P3 


309(3.4*) 


313(7.5) 


P4 


249(2.8*) 


233(5.6) 


P5 


456(5.1) 


64(1.5) 


P6 


251(2.8*) 


2(0.04) 


P7 


219(2.4*) 


24(0.6) 



Table 4: Number of events and the significance for C = 1 fb~ from pure stop-stop production 
using NLO cross sections from PROSPINO. Significance for each case is given in the parentheses 
and entries marked with * (**) indicate that the signal is observable for 1 < C < 5 fb~ (5 < C < 



10 fb 



-l' 



. For this table m~ is 207. 



We have used softer cuts proposed in [14, 20]. The m? mass reach via the above signals 
has been checked using different EW scenarios. 

For the blj fix, I = e and /i, signal we have used the following cuts from [20]. We have 
selected events with one isolated lepton, one tagged b jet, at least 2 jets, fix > 75 and Pt of 
tagged b jet < 80. The SM backgrounds after these cuts are available in Table 6 of [20]. 

For the brj jEx signal we have used the following cuts from [14] which studied this signal 
for the first time. We have demanded events with one tagged b jet, one tagged t jet, no 
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isolated lepton, Et > 70 and Pt of tagged r jet > 40. SM background details after these 
cuts are available in Table 1 of [14]. 

In the light stop scenario if BR (yf 1 ~~ *■ T\V T , u t t)=100 % then it is estimated that rrij < 
280 (305) can be probed by brj Et signal with C = 5(10)fb~ 1 . This is similar to the mj 
reach obtained in [14] where the light EW sector was taken to be as in mSUGRA. Typical 
signal sizes for the benchmark points and the corresponding (S/y/B) are presented in Table 
4 for mj = 207. The sensitivity of the signal to the EW sector is illustrated by this table. 

Models with EW sectors as in P2 and P5 (with large leptonic BR of the chargino) give 
the maximum stop mass reach with blj Et signal. For these models mq < 275 (305) can be 
probed by blj Et signal with C = 5(10)fb~ 1 . Otherwise stop mass reach with C = 5(10)fb _1 
varied between 230-260 (260-290) for the rest of the benchmark points. 

It has been noted that the brj Et signal will be disfavoured if f\ is the NLSP and f\ — Xi 
co-annihilation is the main mechanism for relic density production. Similar degradation of 
the blj Et signal occurs if xi is the NLSP and Xi ~ Xi co-annihilation contributes significantly 
to the relic density production. Both the signal will be degraded if t\ and xi masses are 
close together leading to difficulties in b tagging. 

We next consider the blj /Et and the brj /Et signal in the LSG scenario discussed in 
the introduction following [14]. In this scenario we have treated M3 as a free parameter in 
addition to the soft breaking parameters for the third generation squarks. Here all squarks 
except t\ are assumed to be beyond the reach of the LHC 7 TeV experiments. We use 
different EW sectors in the unconstrained MSSM from Fig. 1. 

The first information required is the lower limit on m-x, if any, in this scenario. As noted 
earlier the limits on m-x in models where all squarks (including t\) are assumed to be heavy 
and gluinos are light are changed significantly in the unconstrained MSSM. The lower limits 
on m-x in PI, P2 and P3 are 450, 500, 620 respectively from jets + Et data. One would 
naively expect the above limits to be valid in the LSG scenario as well since the strong cuts 
of ATLAS eliminates the signal from light stop pairs if any [14] . 

But a more stringent lower bound on m-x in the LSG scenario arises from the blj Et 
signal [3]. Using 1.03 fb _1 data the ATLAS group have excluded gluino masses below 500- 
520 for vnq in the range 125 - 300. The m-x limit is insensitive to mj since the hard ATLAS 
cuts eliminate the events from ti pair production. 

We shall use softer cuts [14, 20] to simulate the above signals. However, we have op- 
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ti 


QCD 


W + lj 


W + 2j 


a(pb) 


85.5 


7.7 x 10 7 


1.4 x 10 4 


5.2 x 10 3 


CI 


28.728 


2.2 x 10 5 


3.1 x 10 3 


9.53 x 10 2 


C2 


14.0519 


1.07 x 10 4 


8.2939 


8.5925 


C3 


3.5662 


1.2465 x 10 2 


1.1499 


1.9352 


C4 


1.8647 


12.9337 


0.5005 


0.8375 


C5 


1.7219 


0.0052 


0.4051 


0.7097 



Table 5: The LO cross-sections (including efficiency) of the SM backgrounds after the cuts(Cl - 
C5) . 

timized the cuts in [20] (see below). Our main aim is to find the stop mass reach for a 
fixed gluino mass consistent with the LHC data. For blj fir signal we have implemented the 
following optimized cuts in succession to enhance the signal to background ratio. 

• We have selected events with one isolated lepton (CI). 

• We have selected events with one tagged h jet (C2). 

• We have demanded events with Pt of leading jet > 120 (C3). 

• Events with missing transverse energy (fir) > 70 are selected (C4). 

• Events with M e ff > 300 are selected (C5). 

The response of the SM backgrounds to these cuts are presented in Table 5. We have also 
studied the brj fix signal, analysed for the first time in [14] using the cuts proposed in that 
paper. 

For m^= 550 the reach in vnq for C = 10 fb~ is displayed in Table 6 for different EW 
scenarios. The highest stop mass can be probed in P2 with our blj fix signal at C = 10 fb _1 
is 375. However, the signal may also be degraded due to the degeneracies discussed in the 
light ti scenario. Tick marks in Table 6 indicate the channel via which the highest mq can 
be probed. It is clear that the signals are sensitive to the EW sector and the two proposed 
signals are complements each other. 
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EW sector 
same as 


m-r 

Si 

reach 


Reach obtained by 


blj fix signal 


brj JS? signal 


PI 


310 


V 


- 


P2 


375 


V 


- 


P3 


330 


- 


V 


P4 


325 


- 


V 


P5 


350 


V 


- 


P6 


335 


V 


- 


P7 


305 


V 


- 



Table 6: Stop mass reach in different EW models for mj = 550 for C = 10 fb in the LSG scenario 
using NLO cross sections. Tick mark indicates the channel for which the mass reach is obtainable. 

3 Conclusions 

The present lower limits on the sparticle masses in mSUGRA obtained by the current exper- 
iments at the LHC are primarily governed by the strongly interacting sparticles (the squarks 
and the gluinos). 

However, strong limits on sparticle masses in the EW sector-much above the correspond- 
ing limits from LEP- emerge due to the correlations among these masses in mSUGRA. For 
example, the mass of the LSP is constrained by m~o ~ 160 [40]. This constraint may have 
important consequences for low mass neutralino DM and the direct search for the neutrali- 
nos [13, 39, 40] and neutralino mass reconstruction [49]. The importance of light neutralino 
scenarios of cosmological importance has been discussed in the context of LHC experiments 
by [50]. 

Since the SUSY breaking mechanism is essentially unknown till this date it is worthwhile 
to revisit the above limits in an unconstrained MSSM with no correlation among the sparticle 
masses in the strong and EW sectors. In order to scan the parameter space we have varied the 
sparticle masses focusing on light EW sectors as described in section 2.1. All soft breaking 
masses in this sector are assumed to be below 200. Some important departures from the 
mSUGRA spectrum which have important consequence for neutralino DM as well as new 
physics search at the LHC are SI) no correlation between the chargino and the LSP mass is 
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assumed, S2) the left sleptons could be lighter than the right sleptons and S3) the sleptons 
belonging to the first two generations could be lighter or heavier than the lighter stau mass 
eigenstate. The present paper is, therefore, an extension of [14, 20] which for simplicity 
assumed the masses of the EW sparticles to be correlated as in mSUGRA. The scenario SI 
can be motivated by mSUGRA type model with non-universal gaugino masses [30] while the 
other two scenarios by models involving non-universal scalar masses [31, 32]. 

We then vary the parameters in the EW sector and delineate the parameter space con- 
sistent with the WMAP data on DM relic density (see Fig. 1). As a consequence of SI) 
neutralinos with mass ~ 35 are found to be consistent with WMAP data as long as the LSP 
is bino dominated or is a mixture of bino and wino. This revives the possibility of neutralino 
annihilation via a nearly on shell Z boson and/or lighter h-boson. Moreover a chargino 
NLSP is a distinct possibility in this scenario and chargino-neutralino co annihilation is a 
viable DM producing mechanism as is readily seen from Fig 1. In S2 with relatively light 
L-type sleptons and sneutrino, neutralino annihilation can proceed efficiently. Moreover, 
LSP - sneutrino co annihilation may become an important DM producing mechanism for 
each LSP mass. This mechanism is especially important for S3 when T\ is relatively heavy. 
None of these mechanisms are viable in mSUGRA in view of the LHC data. 

Using the ATLAS selection criteria [2] we find that in an unconstrained MSSM, the 
limits on the squarks and gluinos derived in mSUGRA changes either way by at most 10- 
12% for most choices of the light EW sector consistent with the WMAP data as described in 
Fig. 1 (sub-section 2.3). This shows that i) the limits on m~ and m^ derived in mSUGRA 
are approximately valid in the unconstrained MSSM and ii) EW sparticles much lighter 
than that in mSUGRA are consistent with the LHC data. If the leptonic BR of chargino 
decays are much larger than that in mSUGRA, then the lepton veto in the ATLAS selection 
procedure relaxes the gluino mass limits appreciably. This is illustrated by the points P2 
and P6 which belongs to the scenario S2. It may be recalled that a similar relaxation was 
noted in the context of old Tevatron experiments (see the 1st paper of [33]). For very heavy 
squarks, the gluino mass limit L4 (see sub-section 2.3) reduces by 20 % . This is the largest 
relaxation in m^ noted by us and is illustrated by the point PI belonging to scenario SI. It 
is, therefore, fair to conclude that the current LHC data in the jets + $t channel is by and 
large insensitive to the choice of the EW sector and EW sparticles having masses just above 
the corresponding LEP limits are very much allowed. 
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In the absence of a suitable e + — e~ collider it not easy to directly test the EW sector 
even if the sparticles in this sector are light. The simulations of the LHC 14 TeV experiments 
[18] suggest that the clean tri-lepton signal from chargino-neutralino production or the di- 
lepton + ^j- signal stemming from light slepton pair production may not be viable at the 
experiments at 7 TeV. 

We have, therefore, considered two scenarios each with a subset of the strongly interacting 
sparticles and a light EW sector within the reach of the LHC, They are i) the light stop 
scenario and ii) the light stop gluino (LSG) scenario. In both cases we consider blj fir [3, 14] 
and brj fir [20] signal. We also follow the selection criteria suggested in these references. 

In the light stop scenario we find that if the lighter chargino decays into r rich final states 
with large BR (see, e.g., points P3 and P4 in Table 2) mj upto 260 (305) can be probed by 
the brj fa signal, investigated for the first time in [14], for £ = 5 (10) fb _1 . On the other 
hand if % decays into eor/i with large BR (see,e.g.,P2 and P5) m-r upto 275 (305) can be 
probed by the blj ^t signal for C = 5 (10) fb~ . For the rest of the benchmark points the 
mass reach varies between 230 - 260 (260 - 290) for C = 5 (10) fb" 1 . It was noted in [14, 20] 
that the signal size is quite sensitive to the EW sector. This is also reflected by Table 6. 
It should, however, be borne in mind that the above signals can be degraded due to some 
unexpected degeneracies in the spectrum as discussed in sub-section 2.3. 

Signal of the LSG scenario in the blj $t channel has also been considered by the ATLAS 
collaboration. However, as pointed out in [14, 20] they use very hard cuts which practically 
eliminate the signal from t\ pair production. As a result the gluino mass reach obtained by 
them is practically insensitive to vnq . Using softer cuts following [14, 20] we find the reach 
in mq for a fixed rn-x= 550 and display the results in Table 6. The signal size indeed show 
some sensitivity to the EW sector. Of course the signal will be degraded if the degeneracies 
as discussed above are present. We conclude by noting that supersymmetry with a light 
electroweak sector consistent with the WMAP and LEP constraints is very much allowed by 
the current LHC data. 
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